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We study a two-dimensional compressible Ising spin glass at constant volume. The spin inter- 
actions are coupled to the distance between neighboring particles in the Edwards- Anderson model 
with ± J interactions. We find that the energy of a given spin configuration is shifted from its incom- 
pressible value, Eq, by an amount quadratic in Eq and proportional to the coupling strength. We 
then construct a simple model expressed only in terms of spin variables that predicts the existence 
of a critical value of the coupling above which the spin-glass transition disappears. 

PACS numbers: 75.10.Nr,75.40.Mg,05.50.+q 



Lattice compressibility affects the nature of phase tran- 
sitions in a variety of spin systems. The 2-dimensional (2- 
D) triangular Ising antiferromagnet, for example, is fully 
frustrated and shows no transition to an ordered state; 
however, when compressibility is added, this system ex- 
hibits a first-order transition to a striped phase 0, 0. 
In the (unfrustrated) Ising ferromagnet, the introduction 
of compressibility changes the transition from second- to 
first-order so that the onset of nonzero net magnetiza- 
tion is simultaneous with a discontinuous change in the 
volume 0, 0] . Frustration is central to the nature of the 
spin-glass transition as it leads to large ground-state de- 
generacies [H E3, 0] ■ Because different states with the 
same spin-glass energy can couple differently to local lat- 
tice deformations, compressibility lifts the ground-state 
degeneracy and may dramatically alter the nature of 
the transition and the low-temperature spin-glass phase. 
Since magneto-elastic effects are always present to some 
degree in physical systems, their inclusion in spin-glass 
models may help explain some of the outstanding puzzles 
in spin-glass experiments We report here the effect of 
compressibility on a 2-D spin glass in which the lattice is 
allowed to distort locally while the entire system is held 
at constant volume. 

The ground states of the previously studied compress- 
ible systems were always states that had already been 
ground states of those same systems without lattice dis- 
tortion. In contrast, we find that as the coupling between 
magnetic interactions and lattice distortions is increased, 
the constant-volume compressible spin glass prefers spin 
configurations which had previously been excited states. 
We further find that the critical region just above the 
spin-glass temperature is suppressed as the coupling to 
lattice distortions increases so that above a certain value 
of the coupling, the spin-glass transition is eliminated 
entirely. 

We study compressible two-dimensional Ising spin 
glasses on square lattices with periodic boundary con- 



ditions at constant volume. The Hamiltonian is 



H = 



JijSiSj+a ^ JijSiSj(rij 

(id) 



-r )+U Xa 



(1) 



The first term is the energy of the standard (incompress- 
ible) Edwards- Anderson spin glass, with a sum over near- 
est neighbors. The spins Si take the values ±1, and 
Jij 6 {±J}. The second term couples the spin inter- 
actions to local lattice distortions; r.y is the distance 
between particles i and j, and r$ is the particle sepa- 
ration of the undistorted lattice. We consider the cou- 
pling to linear order with coupling constant a. This term 
allows particles to move based on their magnetic inter- 
actions with their neighbors: a satisfied interaction (i.e. 
JijSiSj = +1) is enhanced when the particles move closer 
together, and an unsatisfied interaction is mitigated when 
the particles move apart. Because the distortions cannot 
be independent, the energy of a given configuration will 
depend on the local arrangement of satisfied (short) and 
unsatisfied (long) bonds, thus breaking the degeneracy 
inherent to the incompressible model. 

The final term in the Hamiltonian represents the en- 
ergy of the lattice itself: Hooke's-law springs of uniform 
spring constant k connect neighboring particles; springs 
are also required between next-nearest neighbors (across 
the diagonals of the squares) to prevent shear. All springs 
have their unstretched length equal to the natural spac- 
ing of particles on the undistorted lattice. 

We identify two important parameters: 5 = Ja/k is 
a length scale which determines the typical size of the 
distortions from the uncompressed locations; the dimen- 
sionless quantity fx — Ja 2 /k gives the average energy of 
the distortions relative to the original spin-glass energy. 
In our simulations, the parameters are chosen such that 
(5 is held fixed at ro/10 while /i varies over the range of 
interest. The magnetic interaction strength J, which we 
set to unity, sets the overall energy scale. 

We prepare a series of spin states and relax the lattice 
to its minimum potential energy via conjugant-gradient 
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FIG. 1: (a) The energy spectrum of the incompressible spin 
glass consists of a series of delta functions separated in energy 
by 4J. (b)-(d) As fx increases from zero, the spectrum shifts 
downward in energy, and each level spreads into a Gaussian- 
shaped band. The number of states under each curve remains 
constant. These data are from a single L = 20 realization at 
T = 0.65. (e) The average of each band is linear in /x; the lines 
associated with different bands have slightly different slopes. 



FIG. 2: (a) The average energy shift, AE = E - Eo, di- 
vided by fi, plotted as a function of the initial energy; both 
axes have been scaled by the system size, L 2 . Data for dif- 
ferent system sizes approach a common parabolic curve as L 
increases, (b) The scaled data for the width of each energy 
band also approach a common curve for large system sizes, 
but the scaling form indicates that the width is only linear in 
L. In the thermodynamic limit, the spread of each band is 
negligible compared to the average shift in energy. 



minimization [lfj . These states are generated using stan- 
dard Monte Carlo techniques on the incompressible spin- 
glass Hamiltonian. For L = 3, 4, 5, we enumerate all 2 
spin states for each bond realization. For larger systems, 
we acquire data from the ground state up to the highest- 
entropy states (where the spin energy is approximately 
zero) by running at multiple temperatures. We average 
over 100 different bond realizations. 

As a function of the coupling, fi, and tmperature, T, 
we examine the probability, P{E), of finding a state with 
energy E. Typical results are shown in Fig. ^ For the 
uncoupled spin glass with [i = 0, P(E) is a series of 
delta functions separated by 4 J || with heights propor- 
tional to g(E) e~ E l T , where g{E) is the density of states 
(Fig- ma)). As \i increases, each initially degenerate level 
broadens and shifts to lower energy (Fig.njb)-(d)) by an 
amount proportional to fi (Fig.^e)). 

As [i is increased, we calculate the energy shift of each 
configuration, AE = E — Eo, where Eq is the energy 
of that configuration at fj, = 0. For each of the original 
levels (i.e., the 5-functions of Fig. ^a)), we compute its 
average shift, AE, and width, a. Both AE and a are 
quadratic in Eq. The data for AE can be scaled onto a 
common curve by dividing by L 2 , the number of particles 



(Fig. Ha)): 

where A = 0.50 and B = 0.13. The smaller shift for low- 
energy states than for those near E = is to be expected 
since low-energy states, with many more satisfied than 
unsatisfied bonds, cannot distort as effectively as those 
with roughly equal numbers of each. 

The data for the width of each band, a, can be scaled 
onto a common curve by dividing by L (Fig. 0b)): 

where C = 0.31 and D = 0.090. While the energy shift is 
proportional to the area, the broadening is proportional 
only to the linear system size. In the infinite-size limit, 
the spreading of each level will be negligible relative to 
its energy shift. Eq. [21 also implies that the energy gaps 
between neighboring levels vanish when /i sa 4/L, i.e., for 
infinitesimal fi in the infinite-size limit. The previously 
discrete spectrum is thus rendered continuous. 

The empirical form of Eq. and the observation that 
the width of a level can be neglected in the thermody- 
namic limit suggest a simplified model. Because the shift 
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FIG. 3: Simulation results for L = 10 systems using the model 
of Eq. [I] solid lines indicate calculations using the entropy 
expression of Eq. Below v* ~ 0.36, the average bare spin 
energy approaches that of the ground state of the incompress- 
ible system (dotted line) as T is lowered. Above v* , however, 
the zero-temperature value of the average spin energy is given 
by Eq. [3 In this regime, the data may be predicted quite ac- 
curately. 



of each level, AE, is simply proportional to Eq , we re- 
place the last two terms in the original Hamiltonian of 
Eq. ^by a single term that gives that shift explicitly: 



JijSiSj + -jp. 

(i,3) 



Jij SiSj 



(4) 



Here we have absorbed multiplicative factors into the 
coupling constant v ps /i/8. The constant offset has 
been discarded as it does not affect the thermodynamics. 
The first term is again the standard Edwards- Anderson 
spin-glass energy. The second term is the resultant of 
the last two terms of Eq. ^ force balance provides a 
typical distortion that depends upon the total spin en- 
ergy, and this distortion then substitutes into the cou- 
pling term, supplying another factor of the spin energy. 
This Hamiltonian contains only spin degrees of freedom, 
with four-spin, infinite-ranged interactions. The long- 
range nature of this model reflects the network of inter- 
connected springs in the original compressible model. 

The simplified model can be simulated directly and 
studied analytically. A straightforward derivation pre- 
dicts a critical value of v at which a zero-temperature 
transition takes place. From the approximate Hamilto- 
nian of Eq. the energy is E = Eq + j^Eq 2 (note that 
because the energy Eq ~ L 2 the last term does not van- 
ish as L goes to infinity) . The minimum of this function 
with respect to Eg gives the level with lowest total energy 
and depends on v according to 



O.min — 



—L 2 
2v 



(5) 



—L 2 



2_Eo,gr 



ound 



(0) 



the first excited state and the ground state of the incom- 
pressible system have the same energy. For large L, the 
ground state energy per spin, En „n ynr\ /L 2 has been de- 
termined to be approximately —1.4 [llj; thus, v* ~ 0.36. 

Above this value, ground states of the uncoupled sys- 
tem are no longer the lowest-energy states of the cou- 
pled system. This is seen in simulations of the simplified 
model (Fig.|SJ. As T is lowered to zero, the thermally av- 
eraged value of Eq approaches that of the original ground 
state when v < v* but approaches the value predicted by 
Eq. [S] above v* . 

Because they involve only spin degrees of freedom, 
the calculations for this model are standard spin-glass 
Monte Carlo simulations using the energy function given 
in Eq. ^ The presence of v brings neighboring energy 
levels closer together in energy, and the large-^ simula- 
tions equilibrate rapidly. We find that the L-dependence 
of the simulations is somewhat trivial: when scaled by 
the number of particles, computed quantities converge 
quickly to asymptotic values as L increases [l^ |. 

The entropy of the ± J spin glass is well described by 
a quadratic with hyperbolic cosine corrections |l2^ : 



S(Eo) = S - SiEo 2 - S 2 cosh(S 3 E Q ) 



(7) 



For small v, this minimum lies in the non-physical range 



For L = 10, typical values are Si = 2.5 x 10~ 3 , 5 2 = 8.1 x 
10~ 3 , and S3 = 0.056. Using this form, we minimize the 
free energy with respect to E to obtain the average E in 
the thermodynamic limit. The solution can be solved for 
numerically. Predicted values for Eq as a function of T 
are shown as solid lines in Fig. [3 agreement with the data 
from finite-size systems is good for large temperatures 
and values of the coupling above v* . The breakdown of 
the prediction near the ground-state energy results from 
the lack of a low-energy cutoff in Eq. [7\ 

A transition occurs at v* . For larger values of v, states 
of low total energy are also high-entropy states. Thus the 
competition between energy and entropy, on which the 
spin-glass transition depends, vanishes. An order param- 
eter characterizing this new transition is the difference in 
spin energy of the level with lowest total energy and of the 
ground state of the uncoupled system, Bo.min — Eq 

.ground- 

Below v* , this quantity is zero; above, it increases as 
i — \i i.e. linearly in v — v* near the transition. 

The tempcraturc-rcduccd fluctuations of Eq are shown 
versus T in Fig. 0Ja) where the solid lines are calcula- 
tions based on the entropy of Eq. As v increases, the 
maximum value of this quantity increases and shifts to 
lower T, and there appears to be a zero-temperature di- 
vergence that occurs as v — * v*. This divergence in the 
fluctuations of the order parameter occurs due to the en- 
ergy levels moving closer together, becoming equal at the 
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FIG. 4: (a) The temperature-reduced fluctuations in (Eo) 
diverge as T —* when v > v* . Solid lines indicate the curves 
predicted from free-energy calculations, (b) The spin-glass 
susceptibility per spin, as a function of temperature. Dashed 
lines are merely guides to the eye. The sharp rise of xsa 
indicates the onset of critical behavior; the temperature at 
which this occurs is pushed down as v increases until there 
is no longer a critical region. The spin-glass phase is thus 
eliminated. 



critical value v* . The specific heat, which includes both 
energy terms in shows no such behavior and always 

goes to zero as T is lowered. 

The scaling behavior of the spin-glass susceptibil- 
ity has been used to determine the temperature and 
critical exponents associated with the spin-glass transi- 
tion 

dill El- Fig.Efb) shows the spin-glass suscepti- 
bility per spin versus T for various values of v. The in- 
crease of this susceptibility as T is lowered signals the on- 
set of the critical phase that precedes the spin-glass tran- 
sition [J^- As v is increased, this critical phase shrinks, 
implying that the spin-glass transition is being destroyed. 

In the T-v plane, there exists an approximate bound- 
ary between the normal paramagnetic phase and the crit- 
ical phase which signals the approach of spin-glass behav- 
ior. This boundary can be characterized, as a function of 
v, by the temperature at which the temperature-scaled 
order-parameter fluctuations are maximized (see Fig.^J. 
It may also be approximated by the temperature at which 



our predicted value of the average spin energy crosses the 
ground-state spin energy (see Fig- El - In either case the 
phase boundary terminates at the same critical value v* . 

We have limited our study to 2-dimensional systems. 
Preliminary results in 3-D indicate that AE and a also 
depend quadratically on Eq with the same dependencies 
on the system size as were found in 2-D. This would 
suggest a similar result for how the spin-glass transition 
can be destroyed due to compressibility effects in 3-D. It 
is important to note that these results rely on the sys- 
tem being maintained at constant volume. We expect 
constant-pressure systems to behave differently since the 
low-energy levels should distort more effectively than 
higher-energy ones, in contrast to the constant-volume 
system. Finally, since the transition is not dependent 
upon the discrete nature of the energy spectrum, its pres- 
ence is not limited to the ± J model but rather should be 
present in Gaussian models as well. 
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